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Abstract. Multipotent differentiation, where cells adopt one of several cell fates, is a determinate 
and orchestrated procedure that often incorporates stochastic mechanisms in order to diversify cell 
types. How these stochastic phenomena interact to govern cell fate are poorly understood. Nonethe¬ 
less, cell fate decision making procedure is mainly regulated through the activation of differentiation 
waves and associated signaling pathways. In the current work, we focus on the Notch/Delta signaling 
pathway which is not only known to trigger such waves but also is used to achieve the principle of 
lateral inhibition, i.e. a competition for exclusive fates through cross-signaling between neighboring 
cells. Such a process ensures unambiguous stochastic decisions influenced by intrinsic noise sources, 
e.g. as ones found in the regulation of signaling pathways, and extrinsic stochastic fluctuations, attrib¬ 
uted to micro-environmental factors. However, the effect of intrinsic and extrinsic noise on cell fate 
determination is an open problem. Our goal is to elucidate how the induction of extrinsic noise affects 
cell fate specification in a lateral inhibition mechanism. Using a stochastic Cellular Automaton with 
continuous state space, we show that extrinsic noise results in the emergence of steady-state furrow 
patterns of cells in a “frustrated/transient” phenotypic state. 


1. Introduction 

Cell fate determination during developmental processes requires the integration of lineage infor¬ 
mation and signaling cues at specific developmental time points to yield robust, reproducible cell 
fate executions. Most cell fate decisions are determinate and are typically coordinated by waves 
of differentiation. Such “waves of fate” have been observed both during the development of the 
Drosophila visual system 03 as well as in the vertebrate retina 0125]. More specifically, for the 
Drosophila visual system differentiation waves are formed in different parts of the embryo following 
specific morphogenetic furrows in order to establish topographic connections throughout the fly vi¬ 
sual system. Additionally, in retinal development, retinal progenitor cells undergo a series of state 
changes before they adopt a final differentiated fate while a state is defined by the cell competence 
to respond to both intrinsic and extrinsic cues. The synchronization of these multiple differentiation 
waves is mediated by specific signaling pathways. One of the most important signaling pathways is 
the Notch/Delta [6l. ! 10,22.25. 3ChTI]. 

The Notch/Delta pathway represents a juxtacrine signaling transduction mechanism in develop¬ 
mental biology for cell fate decisions, in particular in the nervous system 0E2]. It is also known as 
the “lateral inhibition” or the “lateral specification” mechanism, where the selected cell blocks the 
ability of its neighbours to differentiate. More specifically, in lateral cell fate specification neighbor¬ 
ing cells exchange signals in order to adopt a specific fate, while one “local” winner is selected. In 
this regard, in small domains a single “winner” is produced whereas in larger fields a self-organizing 
“salt-pepper” or chessboard pattern is formed [2]. 

The main factors of the Notch/Delta pathway involve the extracellular domain of the two trans¬ 
membrane ligands Delta and Serrate on the surface of one cell interacting with the extracellular 
domain of Notch receptor at the adjacent cell. The ligand/receptor binding triggers proteolytic 
events that result in the cleavage of the receptor’s intracellular domain. Notch interacts with Delta 
in two modes: activating Delta signaling in neighboring cells (trans-activation) while inhibiting Delta 
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Figure 1. Notch/Delta mechanism: The Notch receptor is activated upon binding 
to the Delta ligand that is anchored to the membrane of an adjacent cell. The combi¬ 
nation of intercellular Notch activation and intracellular Delta deactivation results in 
two mutually exclusive fates of initially equipotent neighboring cells. 


signaling in the same cell (cis-inhibition). Lateral inhibition patterning can emerge when Notch sig¬ 
naling downregnlates Delta levels. Under specific conditions a high level of Notch in one cell will 
drive all of its neighbors to low levels of Notch and vise versa in the Delta case. This generates a 
stable lateral inhibition pattern where each over-expressing Notch cell is surrounded by cells with 
increased Delta concentration (see Figure 1). 

Although most of cell fate decisions are fixed within the developmental process, not all of them 
are determinate. Some cell fate decisions require random selection among alternatives. In agreement 
with determinate decisions, mechanisms acting after stochastic choices can compensate to yield 
robust outcomes HZII3T]. Two stochastic influences can play crucial roles: the external and the 
internal noise. The intrinsic noise is associated with stochasticity involved in the genetic, epigenetic, 
transcriptional or translation regulations. On the other hand, extrinsic noise accounts for all the 
microenvironmental factors such as extracellular matrix or stromal components that can influence 
the cell fate determination process. In this work, we investigate the impact of extrinsic noise on 
the robustness of pattern formation introduced by the Notch/Delta mechanism. Here, the extrinsic 
noise coarse-grains the factors impacting Notch/Delta binding at the cells’ membrane. We show 
that extrinsic noise can result in the formation of spatio-temporal “frustrated furrows” associating 
our results with recently discovered biological evidence. Although the word noise is related with 
negative effects as something undesirable, in biological systems noise represents “randomness”. In 
other words, noise may well be regarded as something desirable. 

The structure of this paper is as follows: In section [2j a thorough investigation of recently de¬ 
veloped theoretical models of lateral inhibition is presented while in section 12.11 a simplified “lateral 
inhibition” model of Collier et al. [1996] is defined. The patterning bifurcation diagram in a two-cell 
population is demonstrated in section 13.11 The fine grained “salt and pepper” pattern formation 
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structure is presented and analyzed in subsection 13.1.11 A stochastic differential equation formu¬ 
lation for Notch-Delta is then derived in section 13.21 In subsection 13.2.11 our simulations results 
highlighting the emergence of a “frustrated” cell state associated with the effect of extrinsic noise 
are analyzed and compared with the deterministic case. In section [3~3l we investigate the preponder¬ 
ance of “frustrated/transient” cell fate mechanism in biological systems. A discussion together with 
suggestions for future work is finally given in section [4j 

2. Modeling Lateral Specification 

Several mathematical models have been developed to analyze lateral specification in different 
biological contexts. The first model on juxtacrine signalling was formulated by Collier et al. [1996] 
considering the role of a ligand, Delta, and its receptor, Notch, and demonstrating that lateral 
inhibition is able to generate fine-grained patterns. There have been other extended formulations of 
the aforementioned model, by Owen and coworkers, either by considering the dynamics of ligands as 
well as free and bound receptors [T9] or by inducing positive feedback realizations for both ligand and 
receptor [TOUTS] , applied to different geometries (strings of cells and square and hexagonal arrays), 
and proving that lateral inhibition can generate patterns with longer wavelength. In a more recent 
paper, O’Dea and King [2011] developed a multi-scale technique in order to construct a continuum 
model for investigating the pattern formation dynamics proposed in Collier et al. [1996]. 

Models dealing mostly with genetic circuits through the nonlinear lateral inhibition mechanism 
have been proposed by Plahte [2001] and Shaya and Sprinzak [2011]. More detailed models of the 
Notch/Delta circuit can be found in Hsu et al. [2006] and Meir at al. [2002], More specifically, Hsu 
et al. [2006] demonstrated that the lateral inhibition mechanism favors cells that are in contact with 
fewer inhibitory signal sending adjacent cells while Meir at al. [2002] investigated the parameter 
regime and how this affect the generation of new patterns through the lateral inhibition mechanism. 

Just recently, several quantitative models have been developed in order to address specific experi¬ 
mental questions associated with the specific pathway |I2][2SlSSj • Koizumi et al. [2012] studied the 
role of lateral inhibitory regulation on cells present at the leading edge (tip-cells) during Drosophila 
trachea development. Both their mathematical and numerical results revealed that Notch/Delta 
mechanism enhances the robustness of the tip-cell selection compared with a system regulated by 
self inhibition. The work by Cohen and coworkers [I2j|35] showed that the development of the rni- 
crochaete bristle pattern on the notum of the fruitfly Drosophila melanogaster is dependent upon 
a long-range interacting population mediated by dynamically extending actin-based filopodia. Cells 
use filopodia-like extensions to gather information from non-neighboring cells. Non-local cell interac¬ 
tions mediated by filopodia dynamics are thought to help cells receive morphogen signals as well as 
collect information regarding the identity of nearby cells that provide survival cues. These filopodia 
dynamics generate a type of structured noise that contributes to the formation of a well-ordered and 
spaced pattern of bristles by the induction of intermittent Notch-Delta signaling. 

An investigation of the role of structural noise (spatial as well as temporal) was presented in a 
more recent work of Cohen et al. HU- They proposed an asynchronous cellular automaton model to 
study a lateral inhibition mechanism in a hexagonal lattice. They showed that in the absence of noise 
they were able to recapitulate the results obtained from previous related continuous models using 
a discrete formulation without the need to invoke diffusible morphogens. Whereas the inclusion of 
noise can lead to stripped or spotted pattern refinement whose thickness is directly dependent on 
the noise signal threshold HU. 

Additionally, Sprinzak et al. BU have studied Notch signaling from a synthetic biology perspective. 
In particular, they highlighted the issue of cis-inhibition effect, namely what happens if Notch and 
Delta bind in the same cell alongside with the trans-inhibition effect. In their seminal work, they 
developed different kinetic models that could reproduce the dynamic behavior of the engineered 
Notch-Delta circuits. Mainly, their models were reduced to three equations in which the dynamics 
between the ligand, the receptor and Notch intracellular domain are investigated. 

It has recently been observed by Formosa and Ibanes [2T] that cis-inhibition enriches cell patterning 
by promoting pattern multistability. More specifically, their work predicted novel characteristics of 
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the Notch-Delta pathway such as the role of ligands that travel through the extracellular space and 
the role of diffusible morphogens on patterns wavefronts. Moreover, they showed the antagonistic 
role of Notch acting either as a lateral inducer or lateral inhibitor in specific developmental processes. 
For instance, they predicted that autoactivation of a specific proneural factor drives Delta upstream 
and is a fundamental component for chicken inner ear robustness 021 

A recent work by Barad and co-workers highlighted that the ability of Delta to inhibit Notch within 
the same cell (cis-interaction) can reduce the sources of errors associated with the lateral inhibition 
process [3J0]. They argue that this beneficial effect of noise complements its role in the generation 
of phenotypic diversity. 

Several elements of the Notch/Delta pathway remain unclear. For example, the diffusible transport 
of ligands through the extracellular space suggest a new way of non-local cell-cell communication. 
Another example is the “bipolar” role of Notch/Delta. Namely, it is known that Notch can have 
both an inductive as well as an inhibitory role in specific organs. A thorough review of the aforemen¬ 
tioned models as well as the latest scientific challenges in pattern formation through the Notch-Delta 
pathway can be found at [20]. 

2.1. Mathematical Model of Lateral Specification. In the current work, we construct a model 
based upon the mathematical framework proposed and analysed by Collier et al. [1996]. The model 
is defined by a set of coupled ordinary differential equations that describe the dynamics of signaling 
mediated Notch-Delta intercellular activation and intracellular inhibition. A pair of ODE’s governs 
the evolution of the levels of Notch and Delta activity in each cell i over time r. The Collier et al. 
[1996] model in dimensionless form is given by: 


iii = f (D) - rii , (la) 

di = g(rii)-di, (lb) 

where D denotes the mean level of Delta activity in the cells surrounding cell i, D = yAj Dj, 

where 91j denotes the neighbourhood of cell i and |91,;| specifies the number of its neighbours. To 
model the lateral inhibition exerted by one cell upon another, Collier et al. [1996] defined the 
functions f,g as: 


/ 0 ) = 
g(x) = 


a + x k ’ 

1 

1 + f3x h 1 


(2a) 

(2b) 


where / is a monotonically increasing function modeling the coupling between adjacent cells and g 
is a monotonically decreasing function representing the inhibitory effect of the bound Notch/Delta 
upon Delta production. The positive parameters a,/3,k,h determine the feedback strength. More 
specifically, the value of the exponent k determines whether the Notch-Delta binding is monovalent 
(k > 1) or cooperative (k > 2). Cooperative reaction means that the Notch receptor has more than 
one binding cites, so the receptor has the ability to bind to more than one ligand molecules while 
in the monovalent case (k = 1) the Notch receptor has only one binding site. Here, we assume the 
monovalent case. In this regard, without loss of generality, we approximate function / = ax, noting 
that this latter formulation coincides with the Taylor expansion of equation (I2aj) for x -C 1 and 
a = a -1 . The reduced model is then given by: 
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Parameter a in equation (j5aj) depicts the normalized strength of interaction between neighboring 
cells. Since no explicit discussion of the biological motivation of the parameter estimation chosen by 
Collier et al. [1996] was presented, we set the parameters: 

h = 5, /3 = 100 (4) 

By using a linear response of Notch production we seek to investigate the immediate effect of Delta 
concentration on the production of Notch activity (in the neighboring cells) by just tuning parameter 
a which served as our bifurcation parameter. Below, we will investigate the dynamics of pattern 
formation. 


3. Results 

In this section, we present the results of our analysis for a deterministic and a stochastic version 
of the model. 

3.1. Bifurcation analysis of the Deterministic case in the absence of External noise. 

Firstly, we determine under which parameter regimes “salt and pepper” patterning is achieved with 
neighboring cells adopting opposing fates. Since the Notch/Delta mechanism under consideration is 
local, analyzing the model equations for two cells with periodic boundary conditions gives insight 
into the period-two patterning behavior. Figure 2 presents the bifurcation diagrams of the Notch 
and Delta activation computed as fixed points of Equations I3al3bl with respect to the parameter a. 
The bifurcation diagrams were extracted by performing an arc-length continuation method [28]. The 
existence of a pitchfork bifurcation for a critical value of a leads to a bistable regime of Notch and 
Delta dynamics associated with a saddle point. 




FIGURE 2. Bifurcation diagram: Bifurcation diagram of the Notch concentration 
with respect to the strength of interaction, a, in a pair of cells. The two stable solution 
curves (solid lines) in the central region correspond to the period-two fine-grained 
pattern, while the unstable state corresponds to the homogeneous state (dotted lines). 

The critical point which initiates pattern formation is a cr = 0.3771. 

It can be seen, that if the strength of interaction is sufficiently small the only stable solution is 
the homogeneous state where all the cells have the same concentration of Notch and Delta. At a 
critical value, a cr , the homogeneous state loses stability and gives rise to a pair of heterogeneous 
steady states where the cells adopt two distinct fates: high-Notch/low-Delta (primary fate) and low- 
Notch/high-Delta (secondary fate). Thus, when a > a cr the two-cells system results in a mutual 
inactivation which creates the lateral inhibition pattern. 
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3.1.1. Numerical simulations. We consider a coupled map lattice model with continuous state de¬ 
fined on a two-dimensional regular lattice L = L\ x L 2 , where L\ = L 2 = {1, ...50} are the lattice 
dimensions, with 2500 nodes (cells) interacting in a von Neumann neighborhood with periodic bound¬ 
ary conditions. The state space of the model is ( rii,di ) 6 [0, l] 2 C 3ft 2 . In order to study the fate 
specification wave we initialize our system by assuming a line of cells with opposite fates (inhibiting 
one another) while the rest of the cells are set to a neutral state (no fate state), in other words having 
zero concentration of Notch and Delta (Figure 3(a)). The aforementioned model can be viewed as 
a Cellular Automaton with continuous state space coupled with synchronous update rules with the 
lateral inhibition process described by equations [Ta] and [Tb] In other words, the transition rules 
associated with lateral inhibition can be represented by the following simple mechanism: the level 
of Notch activation, n, reflects the intensity of the inhibition the cell experiences while the level of 
Delta activity, d, reflects the intensity of the inhibitory signal that each cell delivers to its neighbors. 
More specifically, a cell surrounded by neighboring cells expressing high Notch is assumed to be in an 
inactive state (state 0), while a cell surrounded by neighboring cells expressing low Notch is assumed 
to be in an active state (state 1) I 11 each case we performed simulation for parameter values a > a cr 
aiming to investigate how the fate specification wave propagates in space and time. 

We solve the system of coupled equations EH EH numerically using a Runge-Kutta Method of 4 th 
order for each of the N = 2500 cells. Figure El presents the evolution of the fate wave specification 
when a = 1. As expected, the initial line of fate specific cells evolves through time creating a fate 
wave that results in the evolution of a “salt and pepper” pattern. 




■ Low Notch 
I I High Notch 


(a) 

T = 0 


(b) 

r = 20 




(c) (d) 

r = 30 r = 40 

Figure 3. Deterministic fate specification wave: The generation of a “fine¬ 
grained” chessboard pattern in simulations of the deterministic model defined by the 
system of equations (j3a| and (|3b|) with a — 1. 


I 11 figure 4(a) we investigate the evolution of Notch concentration between two different cells. 
Spatial patterning is generated by a wave front, behind which a regular “salt and pepper” pattern 
forms. Initially, the level of Notch concentration in all cells (except the first line with cells having 
already decided their fate) is chosen close to the stable manifold of the saddle point (see Figure E}. 
Figure 4(b) depicts the phase space of the Notch concentration between two consecutive pairs of cells 
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Figure 4. Impact of differentiation on cell dynamics: (a) Time evolution of 
the Notch concentration of different cells. The vertical lines depict the time of fate 
selection for each cell, (b) Phase space of Notch concentrations between different cells. 
The trajectories correspond to consecutive neighboring couples of cells which adopt 
different fates. The yellow line is the stable manifold of the saddle node (see text for 
further information). 


which finally adopt different fates. The stable manifold of the saddle point (denoted in yellow in 
the plane ( ni,rij ) in Figure 4(b)) separates the attraction of the two different states |38J. The cells 
remain on the stable manifold until the wave fate specification forces them to follow the steady-state 
pattern. 

3.2. Stochastic case: The emergence of the third “frustrated” state as a consequence 
of external noise. Here, we investigate the impact of extrinsic noise on the pattern formation 
dynamics. We assume that extrinsic noise influences the cell’s ability in receiving free Delta D and 
that the parameter a is a random variable a = /i + rj T , with \x e [0,1] and r\ T ~ Af(0, a 2 ). The rest 
of the parameters are defined by equation (J3J) presented in section 2. Thus, the stochastic version of 
the deterministic model defined in equations (l3al l3bj) is: 


rii = ijD — rii + Dr] T 


(I; 


1 + fin? 


di 


(5a) 

(5b) 


3.2.1. Numerical simulations. Following the simulation setting presented in section IT31 we simulate 
the evolution of the fate specification wave by initializing the system as in the deterministic case. 
For the numerical integration of equations (15aft and (15b|) in each cell i we used the Euler-Maruyama 
method with time step dr = 0.01 [33] . 

As expected, the extrinsic noise perturbs the steady state “salt and pepper” pattern found in the 
deterministic case. In this regard, in Figure [5j atr = 20 we note that cells in the middle of the 
domain select a fate before the arrival of the differentiation wave. As time evolves in r = 40, we 
observe that most of the cells have already been differentiated. However, several closed paths of 
“frustrated/transient cells”, i.e. undecided ones, emerge enclosing “chessboard” patterned cells. By 
time r = 100, a combination of a fine grained pattern of “salt and pepper” and “furrows of frustrated 
cells” has emerged. 
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(a) 


r = 10 



(c) 

t = 40 



| Low Notch 
| | High Notch 
Frustrated state 


(b) 


r = 20 



(d) 

r = 100 


FIGURE 5. Emergence of frustrated states: The evolution of emerging “frustrated 
furrows” within the “salt and pepper” pattern in simulations of the stochastic system 
of equations (l5al l5bjh with parameters // = 1 and a = 0.01. 



Figure 6. “Frustrated furrows” dynamics: (a) A close look to a neighborhood 
of a “frustrated furrow” reveals cell fate switching, (b) Phase portrait of Notch con¬ 
centrations between different couples of cells. See text for further details. 

By having a closer look at Figure 6(a), we see that these “frustrated furrows” are formed by cells 
that have adopted a fate which is different from the deterministic case. The exact formation of these 
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furrows is varying for different system realizations. As we can observe in Figure 6(a), as soon as the 
differentiation wave has reached these frustrating furrows a conflict in cell fate determination arises: 
cells in the region within the “frustrated furrows” have opposite order of successive black and white 
nodes compared to cells in the region outside the furrows. Cells located at the vicinity of these two 
regions remain undecided, in terms of their fate specification, and therefore the frustrated/transient 
or mixed state perpetuates in time (frustrated cells denoted by gray in Figure 6(a)). 

Additionally, in Figure 6(b) we highlight the Notch concentration of two neighboring cell pairs. In 
one pair (blue curve) noise forces the trajectory to cross the stable manifold, implying the selection 
of opposite cell fates compared to the deterministic case as shown in Figure 4(b) while the other pair 
of cells (red curve) stay close to the separatrix. Furthermore, in the stochastic model “fate selection” 
takes place earlier than in the deterministic case (Figure 7). 



Figure 7. Effect of noise on cell fate specification: Comparison of a cell’s 
Notch level time evolution in the presence of noise (/x = 1 and a = 0.01) and in the 
deterministic case [a — \x~ 1 and a = 0). Stochastic fluctuations can induce cell fate 
selection before the arrival of the differentiation wave. The vertical line illustrates the 
time of cell fate decision in the deterministic case. 


Finally, in Figure 8(a) and 8(c) we compare the state distribution in the deterministic (fi = 1 
and a = 0) and the stochastic case (p, = 1 and a = 0.01), averaging over 100 simulations for time 
t = 100. The intermediate states denote the existence of frustrated cells found in the “fate switching 
furrows”. Additionally, in Figure 8(b) and 8(d), we compare the Notch dynamics in a cross-section 
of the lattice for both cases. As expected, in the deterministic case half of the cells adopt the high 
Notch fate while the rest adopt the low Notch one. On the contrary, in the stochastic case, the 
emergence of frustrated cell furrows divides the domain into patterned regions. 

We also studied the effect of varying the noise amplitude. A significant range of noise levels 
gave rise to emerging frustrated regions within the “salt and pepper” pattern. By increasing the 
level of noise beyond a > 0.3 we obtained a fully disordered system (Figure [U]). Additionally, we 
investigated the effect of neighborhood topology, i.e. the signaling communication range. For a Moore 
neighborhood (8 immediate neighbors) a refined stripped pattern is obtained (Figure fT0l a)h In a 
similar manner, when noise is implemented frustrated furrows emerge within the stripped pattern, 
as shown in Figure fTOf bh The aforementioned results are in line with those obtained by Cohen et 
al. [2010] using a hexagonal lattice. In the same study, the authors suggested that Notch/Delta 
spatiotemporal dynamics are similar to phenomena found in solid state physics, such as grain growth 
and recrystallization mm- Adding to this discussion point, our results suggest that dynamics of 
the frustrated patterns are similar with spin glasses phenomena. In the next section, we discuss 
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the biological relevance of the Notch/Delta mechanism and the existence of frustrated/transient 
phenotypic cells states. 



FIGURE 8. Notch probability distribution: (a) Notch probability distribution in 
the deterministic case, when a = n — 1. Half of the cells adopt the high Notch fate and 
the other half the low one. (b) Notch level in one-dimensional cross sections for the 
deterministic case, (c) Notch probability distribution in the stochastic case for fi = 1 
and a = 0.01 averaged over 100 simulations for time r = 100. Some cells adopt an 
intermediate state, (d) Notch level in one-dimensional cross sections for the stochastic 
case. 


3.3. The frustrated/mixed state in cell-fate decision making. On the basis of the above 
considerations, we have raised the possibility that during development cell decision making processes 
do not follow a deterministic chain of events in adopting their final fate. Instead, we propose that 
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Figure 9. Disorder for high intensity noise: Simulations of the stochastic system 
of equations fl5al l5b|) . with parameters /i = 1 and a = 0.35. Furrows of cells in a 
“frustrated”/transient state evolve through the system resulting in a fully disordered 
pattern. 



Figure 10. Neighborhood effect on the Notch/Delta model: (a) A pattern of 
consecutive high Notch (low Delta) and low Notch (high Delta) stripes is obtained in 
the case we assume interactions with 8 neighbors for each cell, (b) The emergence of 
frustrated furrows in simulations of the stochastic system of equations (15al l5bl) , with 
parameters /i = 1 and a = 0.02 when each cell interacts with 8 neighbors. 


extrinsic noise in cell fate specification signaling may force cellular sub-population to a metastable, 
transient, frustrated state. Even though many of the components regarding differentiation programs 
are known, how noise may affect cell fate specification over time remains poorly understood. In this 
regard, it is not known if stochastic fluctuations are required before the decision point or after the 
decision has been made. 
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Although a role for Notch/Delta signaling pathway as a lateral inhibitory mechanism is well es¬ 
tablished, it has been a subject to considerable probing of many other of its attributes in retinal 
development. In this regard, apart from its well known role as a cell fate mechanism, Cepko and 
coworkers [251126] proved recently that it could also be recruited to allow a preservation of a pool 
of retinal progenitor cells(RPCs) in a transient/frustrated competent state for later cellular diversi¬ 
fication. During this competent/transient state, cells interact with both inhibitory and stimulatory 
environmental (extrinsic) factors to reach a fate decision. In a similar manner, Shin et al. [H] in¬ 
vestigated the role of Notch signaling in regulating a motoneuron-interneuron fate decision. They 
proposed that transient Notch inactivation disrupts lateral inhibition, committing more precursors 
to lineages that produce ventral spinal cords (also known as pMNs) and interneuron cells. Whereas 
continuous inhibition of Notch activity produces excess pMNs and a deficit of interneurons. 

The notion of mixed cell fate decision echoes recent experimental studies of the cellular mechanisms 
underlying retinal epiphysis. Cau and coworkers m revealed a new functional role for the Notch- 
Delta pathway apart form its seminal role as the binary fate switcher. Namely, they proposed that 
Notch activation is required in order to seggregate epiphysial neurons that are in a transient phase 
of “double” or “mixed” identity. More specifically, in epiphysis there are two main type of cells; 
projection neurons and photoreceptors. However, a new type of cells has been identified of mixed 
identity. These cells have not yet received the postulated photoreceptor inducing signal while they 
have also not yet downregulated the projection neuron signal. Therefore, they remain in a mixed 
identity state. However, these mixed identity cells are rarely observed in wild-type epiphysis. Thus, it 
is speculated that either those cells do not represent a true cell “fate” or due to the lack of appropriate 
cell markers their “mixed” state is difficult to be distinguished. 

Another example of a “mixed” cell phenotype is the intrinsically photosensitive retinal ganglion 
cells (ipRGCs) of the mammalian retina that reflect both the characteristics of projection neurons 
as well as photoreceptors [H[511ISl ; 23j- ipRGCs are able to express the melanopsin photopigment 
while there are also photosensitive. In this regard, ipRGCs represent a potential example of a 
mixed/frustrated cell fate. It is speculated that retinal ganglion cells (RGCs) have lost through 
evolution their photosensitivity characteristic apart from their subfamily of ipRGCs that managed 
to retain this dual characteristic. It would be rather interesting to test biologically whether Notch 
activation can seggregate ipRGCs dual function. 

Our findings draw also parallels to a recent work by Alfonso Martinez Arias and coworkers H6JH7J 
that investigated the intermediate transition states that arise during cell fate decision processes as 
a result of the co-operation of two signaling pathways (Wnt and Notch). This transition state is 
an emerging highly unstable and reversible state between the initial (pluripotent) and differentiated 
state. 

Lately, Vistulo de Abreu and coworkers mm proposed that frustrated cells can provide another 
approach of functional immune system. They considered a system of three cells where they try to 
form stable conjugates. For example, if cells “A” and “B” are conjugated and cell “C” is alone, cell 
“C” can form a new conjugate with either cell “B” or “A”, destroying the previous “AB” conjugate. 
Therefore, “A” becomes a frustrated cell by the presence of cell “C”, and returns to the non-conjugate 
state. Every conjugate cell can be destabilized by interactions with neighboring cells. It will be rather 
interesting to investigate the Notch/Delta effect on the aforementioned immune system configuration. 

Finally, the existence of transient/frustrated cell phenotype has already been introduced in the 
epithelial to mesenchymal transition (EMT) [50] , During EMT, cells gain the ability to migrate and 
invade by loosing epithelial characteristics and gaining mesenchymal attributes. However, between 
the epithelial (E) and mesenchymal (M) states, there exists an intermediate phenotype known as the 
partial (P) EMT state. The P state retains both characteristics of epithelial cells as well as features 
of mesenchymal states. Namely, the P state represents a transient/frustrated phenotype that it 
could be transformed to E or M phenotype respectively. Recent experimental and theoretical models 
revealed that several cell types, starting from the epithelial phenotype, converted to partial EMT 
and then to mesenchymal phenotypes as concentration of transforming growth factor-/? (TGF-/?) was 
increased. 
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The “frustrated states” identified in our model may also be associated with stochastic fate switch¬ 
ing. The effect of noise in stochastic fate switching has already been observed and analyzed in the 
context of genetic circuits. More specifically, stabilizing and destabilizing roles of noise in genetic 
circuits have been reviewed in Eldar and Elowitz [2010] (and references therein). 

The aforementioned models highlight that transient behavior can be regarded as an evolutionary 
novel feature of developmental growth. 


4. Conclusions 

The results and ideas discussed above suggest a general framework of the effect of extrinsic noise 
in the Notch/Delta induced cellular pattern formation. In particular, we have adopted a simplified 
version of the model described in Collier et al. [1996], which can produce a “salt and pepper” cellular 
pattern through the use of a Cellular Automaton formulation. We also investigated the impact of 
the differentiation wave, associated with the Notch/Delta mechanism, on the cell fate specification 
for the stochastic and deterministic case. In the deterministic case, pattern formation is coordinated 
by a wave of differentiation. Whereas on the stochastic case extrinsic noise weakens the effect of the 
differentiation wave. Our analysis showed that even small noise intensities can give rise to “frustrated 
cell furrows”. 

However, the observation of “frustrated states” seems to be more an exception than a rule in biolog¬ 
ical development. Obviously, cells have developed de-noising mechanisms. An example is described 
in Barad et al. [2010] who identified the existence of error minimization mechanisms associated with 
the Notch/Delta pathway. Later, Sprinzak et al. [2010] proposed the cis-trans Notch/Delta pathway 
which includes the exact ingredients for the suppression of stochastic fluctuations. Further design 
ideas for de-noising mechanisms can be expected from the analysis of extensions of the mathematical 
model introduced in this paper. 
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